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Abstract 



We compute the decay rate for the Cabibbo-Kobayashi-Maskawa (CKM)-suppressed electro- 
magnetic penguin decay B — > X d + 7 (and its charge conjugate) in the next-to-leading order 
in QCD, including leading power corrections in 1/m^ and l/jn^ ^ n ^ ne standard model. The 
average branching ratio (B(B — > X d + 7)) of the decay B — > X d + 7 and its charge conjugate 
B -> ~X~ d + 7 is estimated to be in the range 6.0 x 1(T 6 < (B(B -> X d + 7)) < 2.6 x 10" 5 , 
obtained by varying the CKM-Wolfenstein parameters p and r] in the range —0.1 < p < 0.4 
and 0.2 < rj < 0.46 and taking into account other parametric dependence. In the NLL approxi- 
mation and in the stated range of the CKM parameters, we find the ratio Rldj/sj) = (B(B — ► 
X d j))/(B(B — > X s j) to lie in the range 0.017 < R(dj / 'sj) < 0.074. Theoretical uncertainties 
in this ratio are estimated and found to be small. Hence, this ratio is well suited to provide 
independent constraints on the CKM parameters. The CP-asymmetry in the decay rates, de- 
fined as a CP (B - X dl ) = (T(B ^ X dl ) - T(B ^ JT dl ))/(T(B ^ X dl ) + T(B ^ JQy)), is 
found to be in the range (7 — 35)%. Both the decay rates and CP asymmetry are measurable 
in forthcoming experiments at B factories and possibly at HERA-B. 
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1. Introduction 

Electromagnetic penguins were first measured by the CLEO collaboration through the ex- 
clusive decay B — > K* + 7 [|IJ , followed by the measurement of the inclusive decay B — > X s + 7 
0. The present CLEO measurements can be summarized as Q: 

(B(B^X s + -f)) = (2.32 ± 0.57 ± 0.35) x KT 4 , 

(B(B -> AT* + 7)) = (4.2 ±0.8 ±0.6) x 10" 5 . (1) 

Very recently, the inclusive radiative decay has also been reported by the ALEPH collaboration 
with a (preliminary) branching ratio 

(B(H b -> X s + 7)) = (3.29 ± 0.71 ± 0.68) x 10" 4 . (2) 

The quantity (B(B — > X s ± 7)) is the branching ratio averaged over the decays B —>■ X s ± 7 
and its charge conjugate 5 — > X, + 7. The branching ratio in (||]) involves a different weighted 
average of the 5-mesons and baryons produced in Z° decays (hence the symbol Hf, ) than the 
corresponding one given in ([!]), which has been measured in the decay T(45) — > B + B~, B°B°. 

These measurements have stimulated an impressive theoretical activity, directed at improv- 
ing the precision of the decay rates and distributions in the context of the standard model (SM) 
and beyond the standard model, in particular supersymmetry. In the SM-context, the complete 
next-to-leading-logarithmic (NLL) contributions have been painstakingly completed [H] - [14 



and leading power corrections in I /ml |~5], |l^, [l7j and 1 / m 2 c |L8], [L^, |20| have also been calculated 



for the decay rate in B — > X s ±7. This theoretical work allows to calculate the branching ratios 
in the SM with an accuracy of about ±9%, yielding (B(B -> X s + 7)) = (3.50 ± 0.32) x 10" 4 
and (B(Hb — > X s + 7)) = (3.76 ± 0.30) x 10~ 4 , in reasonable agreement with the CLEO and 
(preliminary) ALEPH measurements, respectively. The decay rates in eqs. ([]]) and (0) deter- 
mine the ratio of the Cabibbo-Kobayashi-Maskawa (CKM) |22j matrix elements | V t *Vtb/V c b\. 
Since \V c b\ and \ V t b\ have been directly measured, these measurements can be combined, yielding 



\Vts\ = 0.033 ± 0.007 pl| . The central value of \V ts \ is somewhat lower than the corresponding 
value of \V c b\, \V c b\ = 0.0393 ± 0.0028, but within errors the two matrix elements are found to 
be approximately equal, as expected from the CKM unitarity. 

The interest in measuring the decay rate in B — > Xd ± 7 (and its charge conjugate B — > 
Xd + 7) lies in that it will determine the CKM-Wolfenstein parameters p and rj [^| in a 
theoretically reliable way. Likewise, this decay will enable us to search for new physics which 
may manifest itself through enhanced bd'y and/or bdg effective vertices. These vertices are 
CKM-suppressed in the standard model, but new physics contributions may not follow the CKM 
pattern in flavor-changing-neutral-current transitions and hence new physics effects may become 
more easily discernible in B — > Xd + 7 (and its charge conjugate) than in the corresponding 
CKM-allowed vertices bsj and bsg. Closely related to this is the question of CP- violating 
asymmetry in the decay rates for B — > Xd + 7 and its charge conjugate B — > Xd + 7, which 
may provide us with the first measurements of the so-called direct CP violation in B physics. 
With the weak phase provided dominantly by the CKM matrix elements Vtd and V u b in the 
decay B — > X d + 7, the perturbatively generated strong phases can be calculated by taking 
into account the charm and up quark loops in the electromagnetic penguins, which generate 
the necessary absorptive contributions. This calls for an improved theoretical estimate of 
B(B -> X d + 7) and B(B -»• T d + 7) (hence a CP ) in the SM. 
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In what follows, we shall discuss for the sake of definiteness the decays of the 6-quark 
b — > s + j(+g) and b — > d + j(+g), whose hadronic transcriptions are B — > X s + 7 and 
B — > Xa+'j, respectively, but the discussion applies for the charge conjugate decays B — > X s +7 
and B — ► Xd + 7 as well with obvious changes. When it is essential to differentiate between 
the decay of a B meson and its charge conjugate B, we shall do so. The branching ratio for 
the CKM-suppressed decay B — > X d + 7 was calculated several years ago in partial next-to- 
leading order by two of us (A. A. and C.G.) |2H|| . Subsequent to that, the CP asymmetries 
in the decay rates in the leading logarithmic (LL) approximation were calculated in the SM 
2~5, |2B|, and in some extensions of the SM |27], [2"5 |. Much of the theoretical improvements 



carried out in the context of the decay B — > X s + 7 mentioned above can be taken over for 
the decay B — > Xd + 7. Like the former decay, the NLL-improved and power-corrected decay 
rate for B — > Xd + 7 rests on firmer theoretical ground and consequently has much reduced 
theoretical uncertainty; in particular, the one arising from the scale-dependence which in the 
LL approximation is notoriously large, becomes drastically reduced in the complete NLL order, 
presented here. Hence, the improved theory for the decay rate would allow to extract more 
precisely the CKM parameters from the measured branching ratio B —>■ Xd + 7. Of particular 
theoretical interest is the ratio of the branching ratios, defined as 

in which a good part of theoretical uncertainties (such as from m t , fib, &si etc.) cancel. An- 
ticipating this, the ratio R(dj/sj) was proposed in |24| as a good measure of \V t d\- We now 
calculate this ratio in the NLL accuracy, determine its CKM-parametric dependence precisely 
and estimate theoretical errors. 

The CP-asymmetry in the decay rates, defined as 

«cr(B^X d + 1 ) =nB ^ Xd + i)+T(B ^- + i) (4) 

has not so far been calculated in the NLL precision. We recall that, as opposed to the decay rates 
T(B — > X s +7) and T(B Xd + 7), which receive contributions starting from the lowest order, 
i.e., terms of the form (a"(m&) ln n (myi//m b )), the CP-odd numerator in eq. ([|) is suppressed by 
an extra factor a s , i.e., it starts with terms of the form a s {m b ) (a™ \og n (m w / m b )) . To simplify 
the language in the following, we refer to this statement by saying that the CP-odd numerator 
starts at order a s . This results in a moderate scale dependence of acp, arising from the Wilson 
coefficients which contain a term proportional to a s ln(/i;,/m&) which is not compensated by the 
matrix elements in this order. We show the scale-dependence of acp numerically by varying 
the scale fib in the range 2.5 GeV < fib < 10 GeV. The compensation of this scale dependence 
requires the knowledge of the 0{a 2 s ) contributions in the matrix elements of the operators in 
the Wilson product expansion, which is not yet available. However, it is not unreasonable to 
anticipate that a judicious choice of the scale fib in B decays may reduce the NLL corrections. 
Since the results for the CP-even part, i.e., the denominator in eq. ([|), are known in the LL 
approximation, and with the help of the present work now also in the NLL accuracy, this 
information can be used to guess the optimal scale. We make the choice fib = 2.5 GeV for 
which we show that the NLL corrections in the decay rates become minimal (see Fig. 0). Of 
course, one can not insist that this feature must necessarily also hold for acp- Not having the 
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benefit of the complete 0(ot 2 s ) calculation for a C p, this particular choice of /if, is an educated 
guess based on the inclusive decay rates presented here. 

The branching ratio B(B — > X d + 7) and acp depend on the parameters p and rj and this 
dependence is the principal interest in measuring these quantities. To estimate them, we vary 
these parameters in the range —0.1 < p < 0.4 and 0.2 < rj < 0.46, which are the 95% CL. 



ranges allowed by the present fits of the CKM matrix [29]. In addition to this, there are other 
well-known parametric dependences inherent to the theoretical framework being used here, such 
as a s , m t , nib, mc, a em and B s i. We estimate the resulting uncertainty in the branching ratios 
for B — > X d + 7 and B — > X d + 7 in an analogous way as has been done for BiB — > X s + 7) 
(and its charge conjugate). The resulting average branching ratio {B(B — > X d + ^)) = (B(B — > 
Xd + 7) + BiB — > Xd + 7))/2 and the CP rate asymmetry acp are found to be in the range 
6.0 x 10~ 6 < (B(B X d + 1 )) < 2.6 x 10" 5 , and a CP = (7-35)%, with most of the dispersion 
arising due to the CKM parametric dependence of these quantities. For the central values of 
the fit-parameters |2§ p = 0.11, rj = 0.33 one obtains (B(B -> X d + 7)) = (1.61 ± 0.12) x 10~ 5 
and acp = (11.5 — 16.5)%, where the errors reflect the uncertainties stemming from varying 
Pb £ [2.5 GeV, 10.0 GeV] and the rest of the parameters. The ratio R(d-y/sj), defined in 
eq. (Q), is largely free of parametric uncertainties; the residual theoretical error on this quantity 
is small but correlated with the values of p and r]. Amusingly, the theoretical uncertainty 
on this ratio almost vanishes for the central values of the CKM- fit parameters! (see Fig. ^|). 
However, as the presently allowed CKM-domain is large, one can take the largest uncertainty 
AR^dj/sj) I R(d r y/s r y) = ±7%, which is found for 77 = 0.2 and p = 0.4, as an upper limit on this 
uncertainty. The ratio R{d^/s r y) would then provide theoretically the most robust constraint 
on the CKM parameters. For the ratio itself, varying the CKM parameters in the 95% C.L. 
range, we find 0.017 < R{d^f/s r j) < 0.074, with the central value being 0.046. Hence, even a 
first measurement of this ratio will provide a rather stringent constraint on the (p, if) domain. 
We show this for three assumed values, R(dj/sj) = 0.017, 0.046, 0.074 taking into account 
the theoretical errors (see Fig. |6|). 

Although the shape of the photon energy spectra in B — > X s + 7 and B — > X d + 7 is very 
similar, we think that a measurement of the much rarer decay B — > X d + 7 should become 
feasible at future experiments like CLEO-III and B-factories, because these facilities will allow 
for a good K/w discrimination. 

This paper is organized as follows: In section 2, we discuss the theoretical framework and 
present the salient features of the calculation for the decay rates in B — ► X s + 7 and its charge 
conjugate process and the CP asymmetry in the decay rates. In section 3, we work out the 
corresponding decay rates and CP asymmetry for B — > X d + 7, and the ratio R(d r y/s , y). Section 
4 contains the numerical results and we conclude with a summary in section 5. 

2. Decay rates and CP asymmetry in B — ► X s + 7 and B — > X s + 7 

The appropriate framework to incorporate QCD corrections is that of an effective theory 
obtained by integrating out the heavy degrees of freedom, which in the present context are the 
top quark and bosons. The effective Hamiltonian depends on the underlying theory and 
for the SM one has (keeping operators up to dimension 6), 

AC 8 

n eff (b - S7 (+<?)) = -^-\ t J2c,(p)o t (p), (5) 

V2 i=1 
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where the operator basis and the corresponding Wilson coefficients Cj(/x) can be seen elsewhere 
p| . The symbol X t = V^V^* is the relevant CKM factor and Gp is the Fermi coupling constant. 

The Wilson coefficients at the renormalization scale /if, = 0{m^) are calculated with the help 
of the renormalization group equation whose solution requires the knowledge of the anomalous 
dimension matrix in a given order in a s and the matching conditions, i.e., the Wilson coefficients 
Ci(fjL = mw), calculated in the complete theory to the commensurate order. The anomalous 
dimension matrix in the LL |3U| and the NLL approximation |]TTJ are known. The NLL matching 



conditions have also been worked out in the meanwhile by several groups. Of these, the first 
six corresponding to the four-quark operators have been derived in [BTJ, and the remaining two, 



CV(/i = mw) and C 8 (/i = m w ), were worked out in || and confirmed in [[|, [Tj| and |14 
In addition, the NLL corrections to the matrix elements have also been calculated. Of these, 
the Bremsstrahlung corrections were obtained in |5], [24]] in the truncated basis (involving the 
operators 0\, 2 , and 7 ) and subsequently in the complete operator basis @, 0. The NLL 
virtual corrections were completed in ||10|| . This latter contribution plays a key role in reducing 
the scale-dependence of the LL inclusive decay width. All of these pieces have been combined 
to get the NLL decay width T(B — > X s + 7) and the details are given in the literature |JTT]]- JT^ . 

We recall that the operator basis in 7i e // is in fact larger than what is shown in eq. (|5|) 
in which operators multiplying the small CKM factor X u = V u t,V* s have been neglected. If the 
interest is in calculating the CP asymmetry, then they have to be put back. Doing this, and 
using the unitarity relation A c = —\ t — \ u , the effective Hamiltonian reads 

AG 

n eff (b -> S1 (+g)) = {x t [CMOM + cmom + c^o^) + c 2 ( M )o 2 (/x)] 

-A„ [Ci(/x)(0i„M - Ox(ji)) + C 2 (ji)(OM - 2 (fi))\ + •••}. 

(6) 

In this equation terms proportional to the small Wilson coefficients C3, ...,Cq are dropped as 
indicated by the ellipses. The relevant operators are defined as: 

OM = (s L ^T a c L )(c L YT a b L ) , 0,M = (s L ^T a u L )(u L YT a b L ) , 
2 (yu) = (sl7 a1 c l )(cl7 / "6 l ) , 2u (» = (s L 7 At M L )(w L 7 M 6 L ) , 

g_ 

167T 2 " Xr ' y " **" ±W ° Vr/ 167T 



°M = 7^MvKsLcr, v b R )F^ , OM = ^- 2 m b {^){s L T a a llu b R )G a ^ . (7) 



Note that the Wilson coefficients in eq. (|6|) are exactly the same as those in eq. (|5]). Moreover, 
the matrix elements < S7|0j u |6 > and < s , ~fg\Oi U \b > of the additional operators 0\ u and 2n 
are obtained from those of 0\ and 2 by obvious replacements. 

For our intent and purpose, we write the amplitudes for the processes b — > 57 and b —>■ sjg 
in a form where the dependence on the CKM matrix elements is manifest. The amplitude for 
the first process (including the virtual corrections) can be written as 

AG F 

A(b -» 57) = — -j=- < s^\0 7 \b > tree D(b -> 57) , 

D(b^ S1 ) = \ t (A t R + iA t I )+\ u (A R + ^) . (8) 

It is straightforward to construct the real functions A R , A\, A R and AJ from the expressions for 
the virtual correction in ref. JTU] and the NLL Wilson coefficients in ref. The amplitude 
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of the charge conjugate decay b — > S7 decay is then: 

4G F 



A(b -> 57) 
D(b -f s T ) 



v/2 



< S7|0 7 |6 > 4ree D(6 -> 57) 



A*(^ + af)+A;(A£ + L4?). 
The decay rate for the process b —> sj then reads Q 



(9) 



r(6 -> s 7 ) 

\D(b^ S1 )\ 2 



mlG 2 F a em ^ b ^ s ^ 2 



I A, 



32tt 4 



t \2 



+ |A U | 2 (A U R )' 2 + (A T ) 2 



+ 



2Re(X* t X u ) A R A R + AjA] -2Im(X* t X u 



A R A] 



A\A\ 



(10) 



The order a 2 terms, which are generated when inserting the explicit expressions for the functions 
A R , A\, A R , and AJ, are understood to be discarded. The corresponding expression for the 
b — > 57 decay can be obtained from the preceding equation by changing the sign of the term 
proportional to Im(A£A u ). 

An analogous expression for the decay width V(b — > s^yg) of the Bremsstrahlung process, 
where the CKM dependence is explicit, is also easily obtained from the literature || |4], ||, [7|, [X0f| . 
To get rid of the infrared singularity for E 1 — > 0, we included the virtual photonic correction to 
the process b — » sg, as discussed in these references. Another possibility, which was suggested 
in ref. [IT], is to define the branching ratio in such a way that the photon energy E 1 has to be 
larger than some minimal value E™ m . 

It is customary to express the branching ratio B(B 
semileptonic branching ratio B(B — ► X£i/£), 



X s + 7) in terms of the measured 



B(M+7) = E(^±i) 



B(B^X£u e ). 



(11) 



The expression for T s ; (including radiative corrections) can be seen in refs. |33 



In addition to the perturbative QCD improvements discussed above, also the leading power 
corrections, which start in 1/m 2 , have been calculated to the decay widths appearing in the 



numerator and denominator of eq. ( [iTf ) [|15| , |l6l , 17] . The power corrections in the numerator 
have been obtained assuming that the decay B — > X s + r y is dominated by the magnetic moment 
operator O7. Writing this correction in an obvious notation as 



Y{B_-^ Xs + 7) 
Y\B^X~ S + 1 ) 



1 + 



8 h 



m 



2 ■ 



(12) 



one obtains 5b = l/2Ai — 9/2A2, where Ai and A2 are, respectively, the kinetic energy and 
magnetic moment parameters of the theoretical framework based on heavy quark expansion. 
Using Ai = —0.5 GeV 2 and A 2 = 0.12 GeV 2 , one gets 5b /ml ~ —4%. However, the leading 
order (1/m 2 ) power corrections in the heavy quark expansion proportional to Ai are identical 
in the inclusive decay rates T(B — > X s + , y) and T(B — > Xivg). The corrections proportional to 

(8a s )/(37r), present in ref. pcft , into the term A? R . 



2 Note that we have absorbed the factor F = 1 
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A2 differ only marginally. Thus, including or neglecting the l/m 2 corrections makes a difference 
of only 1% in B(B -> ~X~ S + 7). 

The power corrections proportional to 1/m 2 , resulting from the interference of the operator 



O2 (and 0\) with O7 in B — > X s + 7, have also been worked out [18], [19], |20|. Expressing this 
symbolically as 

one finds 5 c /m 2 c ~ +0.03 0. 

It is convenient to express the branching ratio for I? — > X s +7 in a form where the dependence 
on the CKM matrix factors is manifest: 

|2 |\ 12 



8(B - jr. + 7) = ML A + M D . + ^MM B , + tfflo, . (14) 

"cii \"cb\ \> / cb\ \ v cb\ 



The quantities D a (a = t, w, r, £), which depend on various input parameters such as rn t , m b , m c , Pb 
and a s , are calculated numerically and listed in Table 1. The averaged branching ratio 



(B(B — > X s + 7)) is obtained by discarding the last term on the right hand side of eq. fll4|) . 
The CP-violating rate asymmetry has been defined earlier. In terms of the functions D a it can 
be expressed as: 

Since the function Di in the numerator in eq. flT5| ) only starts at order a s , the complete NLL 
expression for acp requires Di up to and including the 0(a 2 s ) term which is not known. Hence, 
in the LL approximation, a consistent definition of acp is the one in which only the LL result 
for the denominator is retained, i.e., in this approximation one should drop terms proportional 
to D u and D r and keep only the LL result for D t , which is denoted as in the following. 
The expression for acp in this approximation then reduces to 

I A< I L>t 

This is what we shall use in the numerical estimates of acp- Note that d[°\ which is also 
shown in Table 1, is proportional to the square of the LL Wilson coefficient C® ,eS To be 
precise, this Wilson coefficient is obtained from the Wilson coefficients Cj (i = 1, ...,8) at the 
matching scale \l = mw by using the 1-loop expression for a s (/i) in the renormalization group 
evolution. Moreover, notice that the power corrections, which are contained in the functions 
D t and D r in the NLL branching ratio (|I4D , drop out in the LL expression for acp- 

In the Wolfenstein parametrization |[23| , which will be used in the numerical analysis, the 
CKM matrix is determined in terms of the four parameters A, A = sin 8c, p and rj, and one 
can express the quantities A t , X u and \V c b\ 2 in the above equations as [Q (neglecting terms of 
0(A 6 )): 

X u = AX 4 (p-t V ), A t = -AA 2 (l-y + A 2 (p-^)V A c = —X u — Xt, \V cb \ 2 = A 2 A 4 . (17) 
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3. Decay rates and CP asymmetry in B — > X d + 7 and B — > X d + 7 

In complete analogy with the .B — > X s + 7 case discussed earlier, the relevant set of 
dimension-6 operators for the processes b — > dj and b — > d^g can be written as 

H eff (b -> d 7 (+^)) = --^ {6 [C7 7 (A*)Or(^) + C 8 (a*)Ob(a*) + C^O^) + 2 (//)0 2 (//)] 

[dOXO^//) - O^/i)) + C 2 (^)(0 2u (^) - 2 (jm))) + ■••}, 

(18) 

where ^ = Vj& with j = u, c, t. The operators are the same as in eq. (0) up to the obvi- 
ous replacement of the s-quark field by the c/-quark field. Moreover, the matching conditions 
Ciimy/) and the solutions of the RG equations, yielding C» coincide with those needed for 
the process b — > s r y{+g). The power corrections in l/m| and 1/m 2 (besides the CKM factors) 
are also the same for Y(B — > X d + r y) and r(5 — > X s + 7). However, the so-called long-distance 
contributions from the intermediate w-quark in the penguin loops are different in the decays 
B — > X s + 7 and B — > X d + 7. These are suppressed in the decays B — > X s + 7 due to the 
unfavorable CKM matrix elements. In B — > X d + 7, however, there is no CKM-suppression and 
one has to include the long-distance intermediate w-quark contributions. It must be stressed 
that there is no spurious enhancement of the form ln(ra u ///&) in the perturbative contribution 
to the matrix elements < X^IO^B > (i — 1, 2) as shown by the explicit calculation in jlO] and 
also discussed more recently in [34]. In other words, the limit m u —>■ can be taken safely. The 
non-perturbative contribution generated by the n-quark loop can only be modeled at present. 
In this context, we recall that estimates based on the vector meson dominance indicate that 
these contributions are small P5] . Estimates of the long-distance contributions in exclusive 
decays B — > fry and B — > in the Light-Cone QCD sum rule approach put the corresponding 
corrections somewhere around 0(15%) for the charged (B^) decays and much smaller 0(5%) 



for the neutral B decays [36, B7J[. Model estimates based on final state interactions likewise give 



small long-distance contribution for the exclusive radiative B decays [|3~8|| . To take this uncer- 



tainty into account, we add an error proportional to D t in LL approximation, viz. ±0.1 Dt°\ in 
the numerical estimate of the function D r when calculating the branching ratio B(B — > X d + r y) 
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In analogy to eq. (|H ) the branching ratio B(B — > X d + 7) in the SM can be written as 

B(B -> JQ + 7) = A + J^l D u + A + A , (19) 



where the functions A (a = t,u,r,i) are the same as in eq. ([141) . While these functions (or 
some combinations thereof) were obtained in partial NLL approximation some time ago |24j] , 
the complete NLL results are presented here for the first time. For numerical values of these 
functions we refer to Table 1. An expression for the averaged branching ratio (B(B — > X d ±~/)) 
is obtained by dropping the last term on the right hand side in eq. flT9]). 

As the branching ratio (B(B — ► X s + 7)) is very well approximated by (B(B — ► X s + 7)) = 
| A t | 2 A/|Kb| 2 , the ratio i2(rf 7 /s 7 ), defined in eq. (0), can be expressed as follows: 



A t A A t r A At 
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The leading term |£t| 2 /|A t | 2 in eq. (pO[) is obviously independent of any dynamical uncertainties; 
the subleading terms proportional to D u /D t and D r /D t are still uncertain by almost a factor 2, 
but numerically small compared to unity (see Table 1). Also, as we shall see in the next section, 
the allowed values of the CKM parameters provide a further suppression of these terms. Hence, 
the overall uncertainty in R^dj/sj) is small. 

Using again the LL expression for the denominator in eq. @, the CP rate asymmetry can 
be written as 

acp(B^ + 7) = -^S^, (21) 



where stands for the LL expression of D t . In the numerical analysis, we will use the 
following expressions for the quantities £j in eqs. (|i~9D and ( pT|) (neglecting terms of 0(A 7 )): 

£ u = A A 3 (p & = AA 3 (l-p + ^), (22) 

with p = p(l - A 2 /2) and 77 = 77(1 — A 2 /2) @. Note that all three CKM-angle-dependent 
quantities £j start at order A 3 . Inserting these expressions into eq. (EHI), a simple form for the 
CP rate asymmetry is obtained: 



a CP (B^X d + 1 )= '' (23) 



A w [(l - p) 2 + ^ 

4. Numerical Estimates of branching ratios and CP asymmetries 



We now proceed to the numerical analysis of our results. Based on present measurements 
and theoretical estimates, we take the following values for the input parameters: a s (Mz) — 
0.118 ± 0.003, m b = 4.8 ± 0.15 GeV, m c /m b = 0.29 ± 0.02, m t = mt (pole) = (175 ± 6) GeV 
(corresponding to m t {m t ) = (168±6) GeV), B sl = (10.49±0.46)%, = (130.3±2.3). For the 
CKM matrix elements we note that the parameters A and A are rather well determined. The 
parameters p and i] are constrained from unitarity fits. The updated fits, taking into account 
also the lower bound on the mixing- induced mass difference ratio AM S / AM^ > 20.4 yield (at 

±1(7) @, 

A = 0.81 ± 0.057, A = 0.22 

r] = 0.33 ±0.065, p = 0.1 l+g;^, (24) 

where A, being very accurately measured, was fixed to the value shown. Note that the allowed 
range of p is now asymmetric with respect to p = due to the mentioned bound on AM S / AMd, 
which removes large negative-p values. We also note that the recent CKM fits reported in [09[] 
yield an identical range for rj but they find p = 0.156 ± 0.090, which is more restrictive for the 
lower bound on p than the analysis in p9[ , that we use here. 



In Table 1, we give the values of the functions D t , D u , D r and _Dj, evaluated for the central 
values of the parameters a s (mz), m t , nib, a em and the semileptonic branching ratio B s i. The 
other two parameters m c /m(, and pb are varied as indicated. We note that the renormalization 
scale dependence of D t is significantly reduced in the NLL compared to the LL result 
As D u , D r , and Di start at order a s only, their p b dependence is more significant, but their 
contribution to the branching ratio is rather small. 
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/x 6 = 2.5 GeV 


Pb = 5 GeV 


p b = 10 GeV 


m c jnib 




0.131 


0.106 


0.086 


0.27 


n (0) /X 4 


U. 14Z 


n 1 1 /i 

U. 114 




u.zy 


n(0) /\4 

A /A 4 


0.155 


0.125 


0.101 


0.31 


A /A 4 


0.150 


0.147 


0.140 


0.27 


T\ / \ 4 

A/A 4 


0.155 


0.154 


0.147 


0.29 


A/A 4 


0.161 


0.163 


0.157 


0.31 


A /A 4 


0.015 


0.011 


0.009 


0..27 


A/A 4 


0.016 


0.012 


0.009 


0.29 


A /A 4 


0.016 


0.012 


0.009 


0.31 


A/ A 


—0.033 


—0.021 


—0.014 


0.27 


D r /\ 4 


-0.043 


-0.028 


-0.019 


0.29 


D r /X A 


-0.055 


-0.036 


-0.025 


0.31 


A/A 4 


0.056 


0.039 


0.028 


0.27 


A/A 4 


0.062 


0.042 


0.031 


0.29 


A/A 4 


0.068 


0.047 


0.034 


0.31 



Table 1: Values of the NLL functions D t , D u , D r , Di (divided by A 4 ) for the indicated values 
of the scale parameter fi b and the quark mass ratio m c /mb- Also tabulated are the values for 
the LL function D^/W 



For the values of the input parameters given above, the theoretical branching ratio for the 
decay B -> Xl+7 in the SM is calculated by us as B(B -> ^+7) = (3.50±0.32) x 10~ 4 . This is 
to be compared with the recent result in ref. ||40|| , where the central value 3.46 x 10~ 4 is quoted 
for the case in which the factor 1/T s i is - like in the present work - expanded in a s . Taking 
into account, that we use |Aj/V^| 2 = 0.96 in the present CKM framework, whereas in ref. [|Q 



a value 0.95 was used for the same quantity, the results here and in [|(| are in agreement. To 
calculate acp in the decay rates for B — > X s + 7 and its charge conjugate B — > X s + 7, we shall 
use the LL approximation ( |T6"D for acp- For the central values of the parameters we obtain (ne- 
glecting corrections of 0(A 2 )): a C p(B — > X s + 7) = —(2.11)7/%, which gives for 7/ = 0.33 ±0.13 
the following prediction for the decay rate asymmetry: a C p(B — > X s + 7) = —(0.70 ± 0.28)%. 
Note that these numbers correspond to the preferred scale //& = 2.5 GeV. For /x = 5 GeV, the 
asymmetry would be acp(B — > X s + 7) = —(0.59 ± 0.25)%. Thus, the direct CP asymmetry 
in B — > X s + 7 in the standard model turns out to be too small to be measurable. 

We now discuss the decay B — > X& + 7 and the CP conjugated process B — » X& + 7. The 
averaged branching ratio (B(B — > X^ + 7)) strongly depends on the CKM parameters p, rj. 
Taking the central values of the parameters a s (M z ), m b , m c /m b , m t , B s i, a em and /x& = 2.5 
GeV, one obtains the following prediction: 

(B(B ^X d + 7)) = 2.43 [(1 - pf +f} 2 - 0.35(1 - p) + 0.07] x 10" 5 , 

~ 1.61 x 10~ 5 [for (p, 77) = (0.11, 0.33), or (p, ff) = (0.107, 0.322)] . (25) 
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In comparison, the result in the LL approximation for (B(B — > X d + 7)), for the same values 
of the parameters is: (B(B -> X d + 7)) = 1.61 [(1 - p) 2 + r/ 2 ] x 1(T 5 . This gives, (B(B -> 
X d + 7)) = 1.45 x 10" 5 for (p,7?) = (0.11,0.33). The difference between the LL and NLL results 
is ~ 10%, increasing the branching ratio in the NLL case (see Fig. p] for yn> = 2.5 GeV). The 
scale (pb)-dependence of (B(B — > X d + 7)) in the LL and NLL accuracy is shown in Fig. [I], 
fixing all other parameters to their central values. 

In Fig. |2| we give the p dependence of the branching ratio (B(B — > X d + y)) for 77 = 0.20, 0.33 
and 0.46, using p^ = 2.5 GeV and the central values for all other parameters. We note that the 
dependence on 77 is not very marked. The branching ratio is largest for the smallest allowed 
value of p (taken here as p = —0.10) and the largest allowed value of 77 (assumed here as 
77 = 0.46), and may reach a value of 2.6 x 10 -5 . The minimum value of (B(B — > X d + 7)) in 
the SM is estimated as 6.0 x 10~ 6 . 

The ratio R{dy/s r y) in eq. (|20|) can be expressed in terms of the CKM parameters p and fj 
as follows (expanding l/|Ai| 2 in powers of A): 



R(dy/sy) = A 2 [l + A 2 (l-2p)] 

~ 0.046 [for (p, 77) = (0.11, 0.33), or (p, fj) = (0.107, 0.322)] . (26) 



Our prediction for the direct CP asymmetry acp{B — ► Xd + 7), based on the LL result fl2lD 
and for the central values of the input parameters, is: 

0.44?? 

acp(B X d + l){pb = 2-5 GeV) 



(1 - p) 2 + ff ' 

-0.16 [for (p, 77) = (0.11,0.33), or (p, fj) = (0.107,0.322)]. (27) 

The scale dependence of this result is as follows: acp(Hb = 5 GeV) ~ 0.13 and acp(pb = 
10 GeV) ~ 0.12. As argued earlier, we prefer p& = 2.5 GeV to estimate acp, as for this choice 
of the scale the NLL corrections in the decay rates are small. 

In Fig. |3] we show the 77 dependence of the direct CP rate asymmetry for the B — > X d + 7 
decay for p = —0.10, 0.11, 0.25 and 0.40, using again p& = 2.5 GeV and the central values of all 
other parameters. The smallest value of acp(B — * X d + 7) is 7% (for p = —0.10 and r\ = 0.20) 
and may reach as high a value as 35% (for p = 0.4 and 77 = 0.46), as can be seen in Fig. |3|. We 
want to stress that the p-dependence of the branching ratio (B(B — > X d + 7)) and both the p- 
and 77-dependence of acp(B —>■ X d + 7) are very marked. Hence, their measurements will help 
to determine these parameters more precisely. 

To that end, it is important to estimate the theoretical uncertainties in the branching 
ratio, the ratio R(drf/s'y), and direct CP asymmetry in B — > X d + 7 and B — » X d + 7 for 
given values of 77 and p. We estimate these theoretical uncertainties by varying the scale 
Pb G [2.5 GeV, 10.0 GeV] and the input parameters in their respective ±lcr-ranges given earlier. 
The procedure adopted is as follows: Individual errors Aj((£>(£> — > JQ+7)}), Ai(R(dy/sy)) and 
Aj(acp) are estimated by varying each parameter at a time and the resulting errors are then 
added in quadrature, much the same way as it has been done for estimating the theoretical 
uncertainty in the branching ratio BiB — > X s + 7). As mentioned earlier, we add an error 
of ±0.1£) t in the numerical estimate of the function D r in order to take into account long- 
distance effects generated by intermediate w-quarks. The resulting theoretical uncertainty on 
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(B(B — > Xd + 7)) from all the sources is shown in Fig. f| as ±ler bands for the central value 
77 = 0.33 as a function of p. For a given value of p and 77, the theoretical uncertainty is: 
A((B(B ^ X d + i)))/(B(B -> X d + 7)) = ±(6 - 10)% on the branching ratio. This is much 
smaller than the factor 4 dispersion in (B(B — ► Xrf + 7)) due to the p and 77 dependence, shown 
in Fig. |2|. 

The uncertainty in the ratio R(d , y / S7) is even smaller, since the theoretical errors on (£>(-B — ► 
A s + 7)) an d (B(B — > Xrf + 7)) tend to cancel. The residual theoretical uncertainty AR/R is 
correlated with the value of p and 77, which is not difficult to see from the relation in eq. ( p6|) . 
For the central value of the CKM fits 77 = 0.33, this is shown in Fig. |5| where we plot R{d r y/ 'sj) as 
a function of p. Interestingly, the theoretical uncertainties almost vanish for p in the proximity 
of the "best fit" value p = 0.11. The largest theoretical uncertainty R{d n //s n f) in the 95% CL. 
allowed CKM-domain is for the point (p = 0.4, 77 = 0.2) where AR(d r y / S7) / R(dj / sj) = ±7%, 
as the ratio R(d r y/s r y) is smallest there. This study suggests that the impact of the measurement 
of R(dj/s r y) on the CKM parameters will be largely determined by experimental errors. 

It is interesting to see with which theoretical accuracy the Wolfenstein parameters p and 77 
get constrained assuming an ideal measurement of R. To illustrate this, we study the fixed- R 
contours in the (p, 77) plane. Our procedure is as follows: We choose three hypothetical values 
R = 0.017, 0.046, 0.074, emerging from our NLL analysis. For each of these values, we solve 
eq. (|26|) for 77. Fixing all input parameters, except p, leads to a curve in the (p, 77) plane 
(fixed-i? contour). Varying then the input parameters a s (m z ), nib, m c /mb, m t and the scale pb 
(one at a time followed by adding the individual errors in quadrature), leads to a band in the 
(p, 77)-plane for each value of R. In Fig. |6] these bands are shown for the values of R indicated 
above. The unitarity triangle corresponding to the "best fit" solution (p = 0.11, rj = 0.33) 
is also drawn for orientation. One sees again that the theoretical uncertainties are minimal 
(practically vanishing) for the "best fit" solution. 

In Fig. |7|, we show the uncertainty on acp due to the scale variation and due to the input 
parameters as a function of 77 (with fixed p = 0.11). As mentioned, the power corrections drop 
out in the LL approximation. For given values of p and 77, we find: A(acp) / a,cp — ±17%. 
Since the asymmetry acp(B — > X d + 7) itself varies between 7% and 35% (see Fig. ||) in the 
presently allowed range of the parameters p and 77, the residual theoretical uncertainty is not 
a serious hindrance in testing the CKM paradigm for CP violation in these decays. Of course, 
it will be nice to complete the calculation for acp in the NLL approximation, which hopefully 
will reduce the theoretical uncertainty on this quantity considerably. 

5. Summary 

To summarize, we have presented theoretical estimates of the branching ratio (B(B — ► 
Xd + 7)) and the ratio R(d^y / '57) in the NLL approximation, and acp(B —>■ X d + 7) in the LL 
approximation in SM, working out also theoretical errors. Varying the CKM- Wolfenstein pa- 
rameters p and 77 in the range —0.1 < p < 0.4 and 0.2 < 77 < 0.46 and taking into account other 
parametric dependences stated earlier, our numerical results can be summarized as follows: 

6.0 x 10~ 6 < (B(B^X d + -f)) <2.6x 10~ 5 , 
0.017 < R(dj/sj) < 0.074 , 
0.07 < a C p{B ^X d + 7) < 0.35 . (28) 
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The central values of these quantities corresponding to the "best fit" parameters (p — 0.11, rj — 
0.33) are: B(B -> X~ d + 7) = (1.61 ± 0.12) x 10~ 5 , #(^7/57) = 0.046 and a CP {B ^X d + 1 ) = 
(11.5 — 16.5)%, with practically no error on R(d , y/s r y). This ratio is also otherwise found to be 
remarkably stable against variation in the input parameters, with the maximum uncertainty es- 
timated as AR(d^ 1 57) / 'R{drj / 57) = ±7% for (p = 0.4, rj = 0.2). These quantities are expected 
to be measurable at the forthcoming high luminosity B facilities. The CP-violating asymmetry 
acp{B — > X s + 7) in the SM is found to be too small to measure. We emphasize the need to 
complete the NLL- improved calculation for acp(B — > X d + 7). 

We would like to thank Francesca Borzumati and Daniel Wyler for helpful discussion. This 
work was partially supported by INTAS under the Contract INTAS-96-155. 
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Figure 1: Average branching ratio of the processes B — > X d + 7 and B — > X d + 7, plotted as a 
function of the scale pb for the central values of the input parameters nib, rn c /nib, B s i, m t , a em 
and a s {mz)- The solid (dashed) curve shows the NLL (LL) result. 
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Figure 2: The p dependence of the average branching ratio for B — > X d + 7 and B — > X d + 7 
is shown for different values of 77: 77 = 0.46 (dashed curve); rj = 0.33(solid curve); rj = 0.20 
(dash-dotted curve). All three curves correspond to pb = 2.5 GeV and to the central values of 
the input parameters nib, m c /mb, B s i, ni t , a em and ot s {mz). The vertical lines show the ±lcr 
range for p from the CKM fits |29| . 
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Figure 3: r\ dependence of the CP rate asymmetry acp(B — > X d + 7) for different values 
of p: p = —0.1 (dashed curve); p = 0.11 (solid curve); p = 0.25 (dash-dotted curve), p = 0.4 
(long-short dashed curve). All four curves correspond to pb = 2.5 GeV and to the central values 
of the input parameters mj, m c /mb, B s i, m t , a em and a s (mz)- The vertical lines show the ±lcr 
range for 77 from the CKM fits ||29|| . 
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Figure 4: p dependence of the average branching ratio for B — > Xd + 7 and B — > Xd + 7 
for fixed 77 = 0.33. The solid curve corresponds to /x& = 2.5 GeV and the central values of the 
input parameters. The upper and lower dashed curves show the theoretical dispersion due to 
the errors in the input parameters m&, m c /mb, B s i, m t , a s (mz),a em and due to the variation 
of the scale G [2.5 GeV, 10.0 GeV]. The long-distance contribution due to the u-quark loop 
is also included in estimating the errors (see text). The vertical lines show the ±ler range for 



p from the CKM fits [29 




Figure 5: The ratio R{d~f/ 'sj) (in %) as a function of the CKM parameter p for a fixed value 
of n = 0.33. The bands show the theoretical uncertainties following from the error estimates 
discussed in text. The vertical lines show the ±lcr range for p from the CKM fits 
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Figure 6: Fixed- R contours in the (p, n) plane, obtained by varying the input parameters and 
the scale pb G [2-5 GeV, 10.0 GeV]. The three bands shown in the figure correspond to R = 0.017 
(bottom), R = 0.046 (middle) and R = 0.074 (top). The unitarity triangle corresponding to 
the "best fit" solution from the CKM fits [^9] is also shown. 




Figure 7: r\ dependence of the CP rate asymmetry acp(B — > Xa + 7) for fixed p = 0.11. 
The solid curve corresponds to pt = 2.5 GeV and the central values of the input parameters. 
The upper and lower dashed curves show the theoretical dispersion due to the errors in the 
parameters nib, rn c /rrib, B s i, m t , a s (mz), a em and due to the variation of the scale pb G 



[2.5 GeV, 10.0 GeV]. The vertical lines show the ±ler range for 77 from the CKM fits [29 
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